Tunneling of a composite particle: Effects of intrinsic structure 
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We consider simple models of tunneling of an object with intrinsic degrees of freedom. This 
important problem was not extensively studied until now, in spite of numerous applications in 
various areas of physics and astrophysics. We show possibilities of enhancement for the probability 
of tunneling due to the presence of intrinsic degrees of freedom split by weak external fields or by 
polarizability of the slow composite object. 
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Quantum tunneling is a subject of constantly renewed interest, both experimentally and theoretically. The standard 
Oh' textbook approach describes the tunneling process for a point-like particle in an external static potential. Chemical 
and nuclear subbarrier reactions especially in astrophysical conditions, as a rule, involve complex objects with 
their intrinsic degrees of freedom. As stated in Ref. [2j], "Although a number of theoretical works have studied 
tunneling phenomena in various situations, quantum tunneling of a composite particle, in which the particle itself 
has an internal structure, has yet to be clarified." There are experimental data [3j, [J] indicating that at low energies 
the penetration probability for loosely bound systems, such as the deuteron, can noticeably exceed the conventional 
estimates. 

The problems of tunneling and reflection of a composite particle were discussed recently with the help of various 
models [H, B 0> B B E3, El ■ It was stressed that new, usually ignored, effects are important for nuclear fusion and 
fission, nucleosynthesis in stars, molecular processes, transport phenomena in semiconductors and superconductors, 
both in quasi-one-dimensional and three-dimensional systems. The resonant tunneling associated with the intrinsic 
excitation, finite size effects, polarizability of tunneling objects, evanescent modes near the barrier, real fl2l . Il3j 
and virtual [l4| radiation processes are the examples of interesting new physics. Below we consider simple models 
which illustrate how "hidden" degrees of freedom can show up in the process of tunneling leading to a considerable 
' enhancement of the probability of this process. 
£ — , Let the tunneling particle possess two degenerate intrinsic states and the incident wave comes to the barrier in 
£NJ ■ a pure state "up" (it is convenient to use the spin- 1/2 language with respect to the z-representation) . We assume 
one-dimensional motion with the simplest rectangular potential barrier of height Uq located at < x < a. At 
low energy E <C Uo, when the imaginary action kqo, — [2m(Uo — £')] 1 / 2 a is very large, the transmission coefficient 
| To oc exp(— 2noa) is exponentially small. This probability can be exponentially enhanced by a weak "magnetic" field 
i applied in the area of the barrier. We assume that the interaction of this field with the particle is —ha x , where h is 
proportional to the transverse magnetic field. 

Indeed, this field creates the "down" spin component and splits the states inside the barrier according to the value of 
a x . In the z- representation the regions x < and x > a acquire the down component in the reflected and transmitted 
waves, 
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Inside the barrier the two tunneling components have slightly different imaginary momenta, k = [2m(Uo — E) — h] 1 ^ 2 
and k' = [2m(Uo — E) + h] 1 / 2 . Correspondingly, the spinor wave function under the barrier is given by 

/ Ce~ KX + C'e~ K ' x + De KX + D'e K X \ 
Mx) = { Ce - KX _ c , e - K ' x+DeKX _ D , eK ' x ), 0<x<a. (3) 

Performing the matching of the wave function, we find the transmission coefficients: for spin up 
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FIG. 1: Quantum tunneling of a composite particle through a barrier. The barrier might be visible by only one of the particles, 
e.g. the neutron within a deuteron does not see a Coulomb barrier due to presence of another nucleus. 



and for spin down (initially not present) 
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and $' = Assuming small penetrability, na 3> 1, k'gs 3> 1, and ignoring exponentially small terms, we obtain 

T+ = 4k 2 ^— ' (7) 
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If, in addition, the splitting is large enough and e "«e 



Ak 2 n 2 



-2ko, 



(K 2 + fc 2 ) 2 

This should be compared with the transmission without splitting, 
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This is valid if the splitting A = kq — k is large enough in the exponents, 2a A 3> 1, 
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that allows one to neglect the exponents e K a . The validity condition, therefore, is that not only n a ^> 1 but, much 
more strongly, a A 3> 1, which is hard to satisfy with electrons and realistic magnetic fields in the laboratory: 
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for an electron with a = 10 nm, Uq — E = 1 eV, requires K$a « 50, so that it should be \h\/(Uo — E) > 1/25, 
or \h\ > 0.04 eV. However, a similar situation can be realized in the case of a system with the ground state as a 
combination of two configurations slightly split by their coupling; this splitting plays the role of the magnetic field. 

There is no need for an external magnetic field if there exists another intrinsic state of the tunneling object that 
would be able to tunnel with a larger probability. In distinction to the cases of resonant tunneling discussed in the 
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literature the situation is possible when the composite particle energetically cannot be transferred to the state with 
favorable conditions for tunneling. However, even the virtual admixture of such an intermediate state can increase 
penetrability. In such a case outside of the barrier the trace of the evanescent state can exist only in the vicinity 
of the barrier. Similar virtual states emerge with necessity in reflection of a composite particle Q as well as in the 
situation when only one of the constituents interacts with the barrier while the rest of the constituents do not feel 
it [TH, HH, [13] • This happens for example at the Coulomb barrier for a system that contains neutral and charged 
constituents (see figure 1). 

Next we briefly consider the case when the excited intrinsic state with energy E' > E can virtually transfer 
excitation into translational energy (the system then is still under the barrier). The two-component wave function 
with the low component describing the new intrinsic state can be written in a form similar to eqs. (JTJBJ) , with simple 
substitutions 

B i e -ikx ^ B , e \x^ pi e ik(x-a) ^ F , e -X(x-a) > ^ 

where A = [2m(E' — E)] 1 / 2 gives the decay of the virtual wave function in free space. The wave function inside the 
barrier still keeps the form ([3]) with k' now describing the imaginary momentum of the virtual state. 
The transmission coefficient here is found as 
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where, analogously to eq. (O we neglect exponentially small terms. This expression has a simple A-independent limit 
for large excitation energy, A>k, k' , which, under the assumption exp(«/a) <C exp(«;a), gives 
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This limiting result corresponds to the sudden breakup of the incident wave function by the edge of the barrier that is 
equivalent to its instantaneous expansion into two components one of which has a strong enhancement of the tunneling 
probability. 

In a more realistic description, the intrinsic wave function of the slow composite particle will change smoothly along 
its path. We consider a one-dimensional motion of the complex of two particles at positions (x, y) with masses m x 
and m y coupled by their interaction V(x — y) and slowly moving in an external potential, U , 

2 2 

H =^ + ^ + V ( x -y)+ U *(x)+U y (y), (15) 

where we allow the external potential U to act differently on the constituents. Introducing the center-of-mass co- 
ordinate R, relative coordinate r, and corresponding masses M and /i (reduced mass), we come to the stationary 
Schrodinger equation 
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+ V(r) + Ux (R+ r -^f)+ U y (R-^)- ^} r) - 0. (16) 



2M dR 2 2/i dr 2 
At low energy we can use an adiabatic ansatz 

¥(iZ,r)=VW(r;iJ), (17) 

where the internal function <f>(r; R) describing the smooth evolution of relative motion parametrically depends on the 
slow global variable R and satisfies the instantaneous equation 

~W + v(r) + u ^ R+r -¥) + u " ( R - r -Jf)- ^} ^ R > = o- (is) 

A loosely bound state is created by the potential V(r) far away from the barrier. As the motion in the direction of 
the barrier proceeds and the wings of the relative wave function propagate in the region of the potentials U, this wave 
function is gradually evolving along with its energy eigenvalue e(R). It was pointed out in [l(| that the polarizability 
of the tunneling system by the field of the target may increase the penetration probability. 

Multiplying eq. (|16[) by 4>*(r; R), integrating out the intrinsic variable r and introducing new functions 

(4>\4>) — N(R), ((j>\d<t>/dR) =N(R)a(R), (4>\d 2 <p/dR 2 ) = N(R)/3{R), (19) 
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FIG. 2: The effective potential U(R) for the deuteron, as obtained from eq. 1231 The dashed line represents the rectangular 
potential barrier. 



the resulting equation for ip(R) becomes 

„ 2M 
xp + 2a{R)i>' + {3(R)ip + —[E - e(R)]ip = 0, (20) 

nr 

where primes label the derivatives with respect to R. Defining the new function u(R), 

i)(R)=u(R)e-i a ^ dR , (21) 
we reduce the problem to the standard Schrodinger equation, 

u" + ^-[E -U(R)]u = (22) 

with the effective potential 

UR = e(R) -E + ^ [a 2 {R) + a'{R) - (3{R)]. (23) 

Here the energy scale is chosen in such a way that far from the barrier, R — > oo, energy e(i?) coincide with the intrinsic 
binding energy Eq. Note that the solution of eq. [TH]is not normalized. But eq. [T7]is, i.e. J N(R)ip(R) is normalized 
to one. 

The closeness of the continuum level would be dangerous in the form of adiabatic perturbation theory where small 
denominators can arise. But in the form of a differential equation, as formulated above, we do not throw away 
non-adiabatic effects. They are important and they make the wave function to evolve. Of course, real dissociation 
is impossible at low energy because of energy conservation. When the binding energy is small, the particles are still 
correlated and can get together after the barrier. (Even in the continuum their wave function would not be the 
product of two independent plane waves, they are still correlated because of the interaction between them.) Our 
model accounts for these features. 

As an application of this approach, we consider the transmission of a composite particle through a rectangular 
barrier, Fig. 1, a problem discussed for a molecule in the context of condensed matter physics in 0,0- For nuclear 
applications we assume the "deuteron" model, when the rectangular potential of Fig. 1 acts on one constituent 
("proton") only, U x = U p , while its partner, a "neutron", is not influenced by the barrier, U y — 0. The intrinsic 
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FIG. 3: Transmission probability for a deuteron incident on a rectangular barrier (dashed-line) and on an effective potential 
barrier (solid line). Potential and energy parameters are given in the text. 



potential V(r) is also taken as a rectangular well with parameters of depth Vq = 11.4 MeV and width ro = 2 fm 
reproducing the deuteron binding energy Eq = —2.225 MeV. The choice of the barrier parameters, Uq — 5.5 MeV 
and a = ro = 2 fm, implies that at R = the intrinsic well exactly coincides with the barrier and the deuteron is 
practically unbound. The masses used to produce the results in fig. 2 are M = 2tojv and fi = mjv/2, where m^r = 939 
MeV is the free nucleon mass. 

The numerically calculated effective potential (23) is shown in Fig. 2 by a solid line. It is obvious that the barrier 
transmission problem is very different from that for penetration through the original potential U(R), a dashed line. 
At R = the deuteron is barely bound with binding energy e(0) = 0.01 MeV; the non-adiabatic terms [the brackets 
in eq. (|23[) ] are very small at R — 0, so that the height of the effective potential here is only £7(0) « 2.21 MeV, as 
seen in Fig. 2. For a deeper intrinsic potential V(r), the height £7(0) would be closer to the top of the barrier but the 
smearing effect of weak deuteron binding is significant. 

The Schrodinger equation ([22)) with the symmetric effective potential E7(i?) allows, at given energy E, for the 
solutions with definite parity, u±. The solution of the transmission problem given by the incident wave from the left 
is their linear combination. We find u+ and u- by a numerical procedure starting from the center of the barrier, 
R = 0, with 

u+(0) = 1, u' + (0)=l, u-(0) = 0, u'_(0) = l, (24) 

and using an arbitrary normalization of these basic solutions. Then we can compute the dimensionless logarithmic 
derivatives at a remote point R = i? , where U(Rq) is negligible, 

A±=i? — t^-. (25) 

The logarithmic derivatives at the mirror point R = — Rq are — A±. 

Now we can perform the matching at R = ±i?o for the scattering (transmission) problem with energy E — ti 2 k 2 /2M, 
where the wave function is 

e ikR + Be -ikR t -oo< R < -R , 
u{R) = { C+u+(R) + C-U-(R), -Ro < R < Ro, (26) 
Fe lkR , R a <R<oo. 



6 



This gives the reflection and transmission amplitudes, 
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where q — fci?o- The transmission probability is given by 
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Such a calculation contains a continuous transition to global energy exceeding the height of the effective barrier. 
However, then one need to take into account the opening of the breakup channels. 

The results are shown in Fig. 3 in comparison with the simple tunneling calculation ([9|) that did not take into 
account the gradual adjustment of the internal wave function. A considerable enhancement of the tunneling probability 
is evident. 

Of course this is just a simple model that does not pretend to give a realistic quantitative description of tunneling 
for a composite object. However, we believe that it is worthwhile to point out that there exist quantum-mechanical 
effects which are not discussed in textbooks and which could be seen even in such a simplified example. Moreover, 
this simple description allows us to clearly demonstrate physics of the process not overshadowed by cumbersome 
computations. 

We use a trial wave function that corresponds to slow motion of the object as a whole at very low energy when 
dissociation channels are forbidden. In the framework of this variational approach we solve the problem exactly 
taking into account non-adiabatic corrections (derivatives of the function describing slow motion) which are usually 
neglected in a standard Born-Oppenheimer approximation in molecular or solid state physics. This leads to the 
differential Schroedinger-type equation that is solved exactly (see also [HI]). 

In conclusion, we would like to stress that tunneling of composite objects is an important topic, regrettably not 
studied in detail. Numerous applications to nuclear, atomic, molecular and condensed matter physics, as well as to 
astrophysical reactions, make the progress in understanding this problem absolutely necessary. 
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